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Abstract 

We investigate an integrable Hamiltonian modelling a hetero-triatomic-molecular 
Bose-Einstein condensate. This model describes a mixture of two species of atoms 
in different proportions, which can combine to form a triatomic molecule. Beginning 
with a classical analysis, we determine the fixed points of the system. Bifurcations 
of these points separate the parameter space into different regions. Three distinct 
scenarios are found, varying with the atomic population imbalance. This result 
suggests the ground state properties of the quantum model exhibits a sensitivity on 
the atomic population imbalance, which is confirmed by a quantum analysis using 
different approaches, such as the ground-state expectation values, the behaviour of 
the quantum dynamics, the energy gap and the ground state fidelity. 

PACS: 02.30.Ik, 03.65. Sq, 03.75.Nt 
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1 Introduction 



The experimental achievement that led to the Bose-Einstein condensates (BECs), using 
dilute alkali gases at ultracold temperatures [1,2], induced a substantial effort dedicated 
to the understanding of new properties of BECs. In particular, the development of the 
techniques used in the production and manipulation of ultracold atoms and molecules [3] 
has opened the way to a new field, the "chemistry" of ultracold systems, i.e where the 
atomic constituents of the dilute gas may recombine forming molecules. Such molecular 
BEG compounds have been obtained by different techniques [4], for instance, by Feshbach 
resonances [5-7] or photoassociation [8]. There can also occur atom-molecule interactions 
that must be at least three-body in nature [9], bringing up new stimulus and challenges to 
our physical understanding. Experimental evidences for three-body recombinations [10] 
as well as for Efimov states [11] provide a physical ground and stimulus for the search of 
triatomic molecular BECs and for the investigation of their theoretical aspects, which is 
our main interest. 

>From the theoretical point of view, ultracold atomic and molecular systems are 
characterized by their large quantum fluctuations. In this sense, it becomes relevant the 
search for exactly solvable models describing atomic and molecular BEC. Indeed this 
has become a very active field of research [12-20], and the experimental relevance of 
these models is currently a very active research subject [21]. Those solvable models are 
expected to provide a significant impact in this area, a view that has been promoted 
in [22,23]. Increasing evidence and recent results show that multi-atomic systems may 
be interesting and relevant for ultracold atomic-molecular in Bose-Einstein condensates. 
A significant question in this context is whether more complex ultracold molecules could 
be created than simple dimers [24]. Also, due to the more sophisticated nature of the 
control of the interatomic interactions, in the case of triatomic molecules, one expects 
a rich quantum phase structure. Indeed, very recent experimental results confirm the 
existence of heteroatomic bosonic trimers in ultracold mixtures [25] which provide us 
with additional motivation to pursue the present investigation. 

In this paper we analyze an integrable model describing a hetero-triatomic molecular 
Bose-Einstein condensate where atomic BECs can combine (in different proportions) to 
produce a compound with two atoms of the same kind and a third one of a different 
species. Our model, that has been shown to be solvable in [26], includes besides the 
interconversion of atoms to molecules and vice-versa, a linear interaction corresponding 
as the external potential and a bilinear interaction corresponding the scattering between 
atoms-atoms, atoms-molecules and molecules-molecules. We start our analysis of this 
model by a classical treatment where we obtain its phase space determining in particular 
the fixed points. We find that for certain coupling parameters bifurcation of the fixed 
points occurs, and we can determine a parameter space diagram which classifies the 
found fixed points. This diagram is determined for the imbalance of the number of atoms 
which allows us to classify it in three distinct cases. Specifically, when the imbalance is 
equal zero or negative there is a spontaneous appearance of additional boundaries in the 
parameter space (three for the zero case and two for the negative case), some of which 
can be identified with bifurcations of the minimum of the classical Hamiltonian. We also 
perform a quantum analysis, where we study the quantum dynamics and compare with 
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the classical results. Here we are interested in studying the ground state of the model, 
because as actual systems are in ultracold temperature some insight can be obtained 
from the ground state. Furthermore as pointed out before the presence of large quantum 
fluctuations make it interesting to look for the phase structure at zero temperature, the 
quantum phase transitions. In our case we are able to look for signatures of quantum phase 
transition. Here we use two definitions, energy gap and ground state fidelity in order to 
find a quantum phase pre-transition, a term that will be explained later. We observe that 
the critical points are pinning down in completely agreement with the classical analysis. 

The paper is organized as follows. In section 2 we present our integrable model. 
Section 3 is devoted to the classical analysis where the parameter diagram are obtained. 
Section 4 is devoted to the quantum analysis where we show the quantum dynamics and 
a study about quantum phase pre-transition. Section 5 is devoted to our conclusions. 

2 The model 

Let us consider the following Hamiltonian describing a hetero-triatomic-molecular Bose- 
Einstein condensate with two identical species of atoms, denoted by a which can be 
combined to a different type of atom, denoted by 6, to produce a molecule labelled by c. 
In terms of canonical creation and annihilation operators {ajbjCja'' jb'^ ,c^} satisfying the 
usual commutation relations [a, a^] = /, etc., the Hamiltonian reads 

H = Ua^Nl + U,,Nl + U,,Nl + UabNaN, + U^cNaN, + U,,N,N, + (1) 

+A*a^^a + IJ^hNh + IJicNc + Vt{o)a^h^c + c^haa). 

The parameters Uij describe (S-wave scattering, jii are external potentials and Q is the 
amplitude for interconversion of atoms and molecules. iVj are the number operators, i.e 
Na = a^a is the number of atoms type a, iV^ = b^b is the number of atoms type b and 
Nc — c^c is the number of molecules. 

The Hamiltonian acts on the Fock space spanned by the (unnormalized) vectors 

\N,;N,;N,) = {aY^{bY^{cY^\0), (2) 

where |0) is Fock vacuum. 

The Hamiltonian above has two independent conserved quantities 

N = Na + Nt + 3iV„ J^Na- 2iV6, 

where is the total number of atoms and J is the atomic imbalance. It is convenient to 
introduce k = J/N, as the fractional atomic imbalance. Since there are three degrees of 
freedom and three conserved quantities, the model is integrable. More details about the 
integrability of this model, using the Bethe ansatz method, can be found in [26]. In what 
follows we will investigate this model in detail. Below we begin with a classical analysis 
of the model and determine the fixed points of the system. 
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3 Classical analysis 

Let Nj, 6j, j = a, b, c, be quantum variables satisfying the canonical relations 

[6,, Ok] = [Nj, Nk] = 0, [iV,, Ok] = t6,kl. 

We make a change of variables from the operators j, j^, j = a, b, c, to a number-phase 
representation through 

j = exp{iej)^/N'j, j = a, b, c, 

such that the canonical commutation relations are preserved. We perform an additional 
change of variables 

z=^{Na + N,-3N,), 
N 

= -^{'^Oa + Ob — Oc), 


such that z and 9 are canonically conjugate variables; i.e. 

[z, 6] = il. 

In the limit of large we can approximate the (rescaled) Hamiltonian by 

H = ^^t^-^' + 2(« -X)z + (3+{z + c+)^/{z + c.){l-z) COS&], (3) 
where we have defined 

A = AiAUa + Uh + U, + 2Uab-2Uac-Ubc), 

a = A[4(c+ + l)f/„ + (c_ + l)t/fe + (c+ + c_ + 2)f/„fe-(l + c+)t/,e-(l + c_)^ 
3 

+ — (2/ia + /i6-/ic)], 
(3 = A[4UaCl + UbC^_ + Uc + 2UabC+C^ + 2UacC+ + UbcC- 

6 

with 

c_ = 1 - 2/c, c+ = l + k, A = — , 

k = ^, fee [-2,1]. 

We now regard (JS]) as a classical Hamiltonian and investigate the fixed points of the 
system. The first step is to find the Hamilton's equations of motion which yields 

dz dH mN^ , n — — - . fee 

Tt = ^ = -^i- + ^^)Viz + c.)il-z)sm^- 
de dH 4^]A^2 



-[2Xz + 2{a-X) (4) 



dt dz 36 

2{z + c-)(l -z) + {z + c+)(l -z)-{z + c+){z + C-) /6g 

2^{z + c.){l-z) '°AA^ 
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The fixed points of the system are determined by the condition 



Due to periodicity of the solutions, below we restrict to ^ G [0, Nn/S]. It is convenient 
to define the functions: 

f{z) = Xz + a — X, 



, . _ 2(1 - z){z + C-) + (1 - z){z + c+) -{z + c+){z + C-) 

~ V(l-^)(^ + C-) 

Note that the domain of g{z) is z G [-1, 1) if /c G [-2, 0] and ^ G (2A; - 1, 1) if A; G (0, 1). 

We observe that the fractional atomic imbalance k plays an important role in the 
behaviour of the g{z) function. For < 0, g{z) is divergent only at 2; = 1, while for the 
case of > 0, g{z) is divergent at 2; = 2A; — 1 and z = 1. Since k affects the domain and 
the shape of the function g{z)^ this property will affect the type of solutions of (I5|). In 
Fig. [Uwe illustrate the behaviour of the function g[z) for different values of k. It is, in 
fact, necessary to treat the cases of A; < 0, A; = and A; > separately. 




Figure 1: (Color online) The behaviour of the function g{z) for three different values of 
k. 



3.1 Negative case : — 2 < A; < 

Here the domain of ^'(2;) is z G [—1, 1) and g{z) is divergent at 2; = 1, but finite at z = — 1. 
This leads to the following classification for the solutions of ([5]): 

• ^ = and 2; is a solution of 
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which can admit zero, one or two 
6 = Nn/6 and z is a solution of 



solutions. 



(6) 



/(^) = -9iz), (7) 
which can admit zero, one or two solutions. 

z = — c+, which vanishes the first equation of ([4]) and reduces the second eq. of ([41) 
to the expression 

^=kT-2^ 2(fc + 2) '''^N^^ 
such that 6* is a solution of 

66 2v/-3fc(fc + 2) a \ 

'''^N^ = 3k - ^ j ' 



for which there are two solutions for I ^3k^^ ~ (-^ ~ ■fc+2) I ^ 1- 

3.2 Zero case: k = Q 

Now we consider the case k = 0, where the domain of g{z) is ^ G (—1, 1) and g{z) is 
divergent at z = 1, but finite at z = —1, similar to the previous case. This leads to the 
following classification for the general problem: 

• = and z is a solution of 

f{^)=9{^), (10) 
which can admit zero, one or two solutions. 

• 6 = Ntt/6 and z is a solution of 



/(^) = -9{z), (11) 
which can admit zero, one or two solutions. 

• z = —1, which vanishes the first eq. of ^ and reduces the second eq. of ([4]) to the 
following linear equation between the coupling parameters 

A =2' (12) 

which can admit just one solution. This result is compatible with that obtained in 
the previous case by taking the limit — in eq. (JH]) 
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3.3 Positive case: < A; < 1 



In this case the domain of g{z) is z E {2k — 1, 1) and g{z) is divergent at both extremes 
of the interval, z = 2k — 1 and z = 1. Now, a different scenario emerges, compared to the 
previous two cases. This leads to the following classification for the general problem: 

• ^ = and z is a solution of 

f{^)=9{^). (13) 
which can admit one, two or three solutions. 

• 6 = Nn/6 and 2; is a solution of 

/(^) = -9i^), (14) 
which can admit one, two or three solutions. 

We can collect all different types of solutions of eq. ([5]) in a parameter diagram, 
dividing the parameter space into different regions, for each case of k discussed above. 
For example, for the case of k positive, to construct this diagram, we observe that the 
boundaries between each regions occur when / is the tangent line to ±g; i. e. for values 
of A and a such that 



for some zq. This requirement determines the boundaries in the parameter space, which 
are depicted in Fig. [2](c) for k = 0.5. 



As in the case of k positive, we can determine the region boundaries in the parameter 
space for the other two cases. However, because of the existence of solutions of the form 
given by (jHl), which do not have an analogue for positive k, we see the appearance of 
new boundaries given by the conditions A = (a =F 9{—k — l))/(^ + 2) for negative k and 
A = a/2forA; = 0. The boundaries in parameter space are illustrated in Fig. [2](a) and 
Fig. [2](b) for k = —1 and k = 0, respectively. Notice that the two additional boundaries, 
which delimit region C, for k = —1 are reduced to a unique boundary for A; = 0, which is 
not present for k = 0.5. Therefore, we have a different scenario for the parameter space 
diagram, depending if the fractional atomic imbalance k is negative, zero or positive, as 
illustrated in Fig. El Basically, we can summarize the typical behaviour of the parameter 
space diagram as follows: when k is negative, the parameter diagram is divided in five 
regions: in region A there is no solution for z when 6 = and one solution for z when 
6 = Nit/6. In region B there are two solutions for z when ^ = and one solution for 
z when 6 = Ntt/G. In region C there is one solution for z when ^ = 0, one solution for 
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a a a 



(a) (b) (c) 

Figure 2: (Color online) Parameter space diagram identifying the different types of solu- 
tion for eq.(l5l) for different values of k = — 1;0;0.5. We observe (a) five distinct regions 
for the negative case; (b) four distinct regions for k = 0; (c) three distinct regions for the 
positive case. In (a) the boundaries are given by A = (a =F g{—k — i))/{k + 2) while in 
(b) it is given by A = a/2. 

z when 6 = Ntt/6 and two solutions for 6 when z = —k — 1. In region D there is one 
solution for z when ^ = and two solutions for z when 6 = Ntt/6. In region E there 
is one solution for z when 9 = and no solution for z when 9 = Nn/G. For the case 
k = 0, region C disappears and the phase diagram is left with the four regions A, B, D, 
E discussed before. When k is positive the diagram is divided in three regions: in region 

I there is one solution for z when ^ = and one solution for z when 9 = Ntt/G. In region 

II there are three solutions for z when ^ = and one solution for z when 9 = Ntt/G. In 
region III there is one solution for z when 9 = and three solutions for z when 9 = Nir/G. 
It is interesting to mention that the fractional atomic imbalance also plays an important 
role in hetero-diatomic molecular Bose-Einstein condensates [27-29]. 

To help visualize the classical dynamics, it is useful to plot the level curves of the 
Hamiltonian ([3]). Since the fixed points change the topology of the level curves, qualitative 
differences can be observed between the different regions. The results are depicted in Fig. 
[3] for A; = —1 (on the left), = (in the middle) and k = 0.5 (on the right). For clarity, 
we use convenient intervals for 9 and z. 

In Fig. [3](a) we show the level curves of the Hamiltonian ([3]) for /c = — 1, illustrating 
the typical behaviour for regions A, B, C and E (from the top to the bottom). In region A 
there are local minima at 69/N = ±7r. Besides the minima at 69/N = ±n, two additional 
fixed points (a maximum and a saddle point) are apparent in region B occurring at 6' = 0. 
In region C there are minima at 69/N = ±7r and for ^ = just one fixed point, a 
maximum. There are also saddle points when z = 0. In region E just one fixed point, a 
maximum, occurs for ^ = 0. 

In Fig. [3](b) we show the level curves for = for the same regions illustrated in the 
previous case, except that now instead of region C there is just one straight line separating 
regions B and D. The behaviour here is analogous to the previous case of negative fc, with 
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the emergence of a maximum (minimum) when passing from region A to B (E to C). 

In Fig. [3](c) we present the level curves of the Hamiltonian ([3]) for A; = 0.5, illustrating 
the behaviour of regions I, II, III and I (from the top to the bottom). In region I there is 
a maximal point at = and a minima at Q6/N = ±n. Two additional fixed points, a 
saddle and a maximum occur, in region II at ^ = 0, while two additional fixed points, a 
saddle and a minimum, occur in region III at 66 /N = ±7r compared to region I. 

We observe that the pattern of the level curves is distinct for the cases of k negative 
and zero compared to the positive case. 

In the following sections we will conduct an analysis of the quantum Hamiltonian. 
We will focus our attention on the case A = 0, in this way the model has one effective 
parameter, a. In particular we will establish that the bifurcation occurring at (a, A) = 
{—g{—k — 1),0) for the negative case and (a, A) = (0.5,0) for /c = can be seen to 
infiuence the ground state properties of the quantum system. 

4 Quantum analysis 

We now turn our attention to a quantum treatment of the model, to investigate the nature 
of the additional threshold couplings for the cases where the fractional atomic imbalance 
k is negative and zero. In particular we analyze the Hamiltonian in the no scattering limit 
where Uij = for all i,j = a, 6, c, 

H = HaNa + fibNb + iJcNc + fi(a^a"^6^c + c^baa). (15) 

This simplifies substantially the Hamiltonian, however it remains sufficiently non triv- 
ial to enable us to gain an understanding of the quantum behaviour through the quantum 
dynamics, ground-state expectation value, gap and fidelity. The no scattering limit corre- 
sponds to the coupling A = in the classical analysis of section 3. With reference to Fig. 
[2] there are two threshold couplings when k is negative and three threshold couplings for 
k = 0. For the case of k negative, one occurs at (a, A) = {—g{—k — 1), 0), signifying the 
bifurcation of the global minimum of the Hamiltonian, while the other occurs at (a, A) = 
{g{—k — 1), 0), signifying the bifurcation of the global maximum. For the specific example 
of A; = —1, these thresholds are (a, A) = (0.866,0) and (a, A) = (—0.866,0), respectively. 
For the case k = 0, there are three bifurcations at (a, A) = (—0.5,0), (0,0), (0.5,0). The 
case (a, A) = (—0.5, 0) signifying the bifurcation of the global maximum, (a, A) = (0.5, 0) 
signifying the bifurcation of the global minimum while {a, A) = (0, 0) signifying the bifur- 
cation of the saddle point. In contrast, there are no bifurcations along the line A = for 
the positive case. We focus our attention to the coupling (a, A) = (0.866,0) for k = —1 
and (a, A) = (0.5, 0) for /c = 0, as in these cases the bifurcation of the fixed point in phase 
space is associated with the ground state of the quantum system. 

4.1 Quantum dynamics 

In general the time evolution of any state is given by |^E'(t)) = U{t)\(f)), where U{t) is the 
temporal operator U{t) = ^^^g |m)(m| exp(— ii^m^), \^) is an eigenstate with energy 
Em and |0) represents the initial state with N = Na + Nb + 3Nc. We adopt the method of 



9 



directly diagonalizing the Hamiltonian (fTSl) as done in [31,32] and compute the expectation 
value of z{t) through. 

im = ^(^WliVa + m- 3iV,|v]/(t)). (16) 

In our analysis, for fixed total number of atoms N and fixed atomic imbalance J, 
we will use the initial state configuration |0, — J/2, {2N + J)/6) for the cases where k is 
negative and zero and |J, 0, (A^ — J)/3) for the case where k is positive. We therefore 
compare the three cases of the quantum dynamics, with fractional atomic imbalance 
negative, zero and positive. 

Results of the expectation value for z are shown in Fig. [4] for the cases of /c = — 1, 
and 0.5. We are using N = 900 and J = -900; 0; 450 for k = -1; 0; 0.5, respectively. We 
fix the parameter = 1 and use /ic as the variable coupling parameter. In terms of the 
classical variables, this corresponds to vary the parameter a. The qualitative differences 
are quite apparent. In the case of k = —1, Fig. [4](a), we find that for a < 0.866 there are 
irregular oscillations in z. Similar behaviour occurs for a < 0.5 for k = 0, Fig. [4](b). As 
the coupling parameter a is increased across the threshold value at a = 0.866, for k = —1 
and a = 0.5 for k = 0, the transition to localized oscillations is significant in cases (a) 
and (b). By comparison the dynamics in Fig. [4](c) for k = 0.5 show a collapse and revival 
of oscillations. 



4.2 Ground state expectation values 

Now using the equation (fTHll , we compute the normalized ground-state expectation value 
3{Nc)/N for the quantum system as the coupling is varied. Results are shown in Fig. [5l 

In general, agreement with the classical result is found: As the threshold coupling 
a = 0.866 (for k = —1) and a = 0.5 (for /c = 0) is crossed, the maximal possible number 
of molecules that can be formed for each case (100% for k = —1 and 50% for A; = 0) is 
reached. In both cases, there is an abrupt change in the expectation value 3{Nc)/N at the 
threshold point. However, for k = 0, the expectation value 3{Nc)/N does not exhibit any 
sudden change, indicative of the fact that there is no boundary in Fig. EJc). Therefore, 
qualitative changes are observed between the cases of k negative and zero and the case of 
k positive. 

4.3 Quantum phase transitions 

In order to gain a better insight into the effect of the threshold couplings for the quantum 
system, in our final analysis we investigate the existence of quantum phase transitions in 
our model (fTSll . 

A Quantum Phase Transition (QPT) is usually defined as a phase transition in the 
ground-state of the system under the variation of some parameter. Basically, there is 
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a sudden change in the structure of the ground state at the QPT, and the properties 
such as entanglement, correlations, etc reflect this sudden change [30]. There are difl'erent 
methods to determine a QPT. In particular, we will study the behaviour of the energy gap 
and fidelity of the system to identify a QPT. Here we mention that a QPT is rigorously 
defined in the thermodynamical limit N ^ oo. For large but finite the system does 
display an increasing sharp distinction between ground state regions, called Quantum Pre- 
Phase Transitions (QPPT). The occurrence of a QPPT in a finite system is a precursor 
for a QPT in the thermodynamic limit. Let us now study the QPPT of the Hamiltonian 

4.3.1 Energy gap 

One possibility to identify a QPPT is through the energy gap, which is defined as the 
difference between the first excited state and the ground-state of the system. 

/\E = El — Eq. 

In Fig. [6] we plot the gap against the coupling parameter a for the cases of /c = —1,0 
and 0.5 using f2 = 1 and different values of N. In all cases the energy gap exhibits a 
minimum, which is much more pronounced in the cases of k negative and zero compared 
to the case where k is positive. We observe that as long as N increases, the point where 
the gap tends to vanish corresponds to a = 0.866 for k = —1 (Fig. [6](a)) and a = 0.5 for 
k = (Fig. [6](b)), in agreement with the classical analysis. In contrast, when k is positive 
there is no abrupt variation of the energy gap as shown in Fig. [6](c) and QPPT are not 
expected. 



4.3.2 Fidelity 

Another possibility to investigate the QPPT is through the behaviour of the fidelity, 
which is a concept widely used in the Quantum Information Theory [33-35]. The fidelity 
is basically defined as the modulus of the wavefunction overlap between two quantum 
states. Assuming the ground state of the system is non-degenerate, let "if (a) denote the 
unique normalized ground state. For fixed small A we define the function fidelity or 
ground-state wavefunction overlap Fid/^{a) by 

Fid^{a) = |(vl/(a(l - A))|xl/(a(l + A)))|, 

which is symmetric in A, bounded between and 1 and satisfies Fido^a) = 1. For systems 
which exhibit a quantum phase transition in the thermodynamic limit, the wavefunction 
overlaps between states in different phases go to zero in this limit. The occurrence of a 
minimum in the ground-state wavefunction overlap in a finite system is then a precursor 
for a quantum phase transition in the thermodynamic limit. Thus for finite systems we 
identify quantum phase pre-transitions at a coupling a for which the fidelity is (locally) 
minimal. Fig [7] shows the behaviour of the fidelity for (a) k = —1; (b) = 0; (c) = 0.5 
for fixed A^(A) and different values of A(A^) on the top (bottom). It is clear that the 
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minimum value of Fid/\{a)^ which determines the quantum phase pre-transition, tends 
to occur at a ~ 0.86 for k = —1 and a ~ 0.5 for A; = 0. The distinction between 
the predicted threshold coupling and the observed pre-transition coupling is that the 
pre-transition coupling also occurs for k positive, although for fixed N the minimum of 
Fid^{a) is substantially more pronounced for k negative and zero compared to k positive. 
In all instances the value of minimum decreases with increasing N. We remark that the 
value of a at which the minimum occurs is independent of A, as shown in FigjT] (on the 
top). In our previous classical analysis qualitative differences are only found precisely 
when k is negative or zero. Here it is clear that the distinguishability of two phases is 
more reliable also for k negative or zero. We then interpret these results as the emergence 
of quantum phase boundaries for k negative or zero. 



5 Conclusion 

We have considered a model describing a mixture of two species of atoms in different 
proportions which can combine to form a bound molecular state at zero temperature. 
This hetero-triatomic molecular Bose-Einstein condensate model has been investigated in 
detail through a classical and a quantum analysis. 

We have found that the fractional atomic population imbalance k, an extra control 
"knob" characteristic to heteronuclear models, plays an important role in the determi- 
nation of the phase boundaries in the diagram of parameters in the classical analysis. 
This property also holds at a quantum level by inspecting the ground-state expectation 
values and the character of the quantum dynamics of the model. We have also looked 
for the quantum phase pre-transitions in our system and shown that the quantities en- 
ergy gap and ground state fidelity are suited for revealing QPPT and pinning down the 
critical (bifurcation) points. 
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(a) k = -l (b) k = (c) k = 0.5 
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Figure 3: Level curves for the Hamiltonian ([3]), where the dark regions indi- 
cate lower values than the light regions. Here we are using for: (a) k = —1 
on the left (A, a) = (0; -1.0), (5; 2.5), (0, 0) and (0,1.5); (b) A; = in the mid- 
dle (A, a) = (0,-l),(2.5,2.5),(0;0) and (0,1.5);(c) k = 0.5 on the right (A, a) = 
(0; -1.5), (5, 2.5), (-5, -2.5) and (0, 1.5). 
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Figure 4: Time evolution of the expectation value of z for the Hamiltonian ( ITSll with 
= 900, for (a) k = —1 and initial state |0,450, 150). We are using natural units. The 
oscillations are largely irregular with significantly decreasing amplitude as the point at a = 
0.866 is crossed. This point corresponds to the boundary at (a, A) = (0.866, 0) between 
regions C and E as shown in Fig. [2](a); (b) A; = and initial state |0, 0,300). A similar 
behaviour occurs as the point at a = 0.5 is crossed. This point correspond to the boundary 
at (a. A) = (0.5,0) as shown in Fig[2](b); (c) k = 0.5 with initial state |450, 0,150). The 
oscillations display collapse and revival behaviour with smoothly decreasing amplitude. 
Here there is no abrupt behaviour, indicative of the fact there is no boundary at A = in 

Fig. EKc). 
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Figure 5: (Color online) Normalized ground-state expectation value of the molecular 
number operator (Nc) versus the coupling parameter a for the three different cases k = —1, 
and 0, 5. Here we are using = 1 and = 900. For the cases k = —1 and k = there 
is an abrupt change in the expectation value 3{Nc)/N as the threshold coupling a = 0.866 
(for k = —1) and a = 0.5 (for A; = 0) is reached. In contrast, for k = 0.5, the expectation 
value 3{Nc)/N increases smoothly with a, not exhibiting any abrupt behaviour. 
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Figure 6: (Color online) Energy gap between the first excited state and the ground state 
as a function of a for (a)fc = —1; {h)k = 0; (c)fc = 0.5 and different values of A^. We are 
using n = 1 and natural units. 
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Figure 7: (Color online) Ground-state wavefunction overlaps as a function of the coupling 
parameter a for (a)fc = —1; {h)k = 0; {c)k = 0.5 and f2 = 1. On the top we are using 
= 900 and different values of A. In the bottom we are using A = 0.01 and different 
values of N. In all cases the fidelity exhibits a minimum, which is substantially more 
pronounced for k = —1 and k = 0, compared to k = 0.5. 
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